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Abstract 

The Fermi-Pasta-Ulam (FPU) lattice with periodic boundary conditions and n particles 
admits a large group of discrete symmetries. The fixed point sets of these symmetries 
naturally form invariant symplectic manifolds that are investigated in this short note. 
For each k dividing n we find k degree of freedom invariant manifolds. They repre- 
sent short wavelength solutions composed of k Fourier-modes and can be interpreted 
as embedded lattices with periodic boundary conditions and only k particles. Inside 
these invariant manifolds other invariant structures and exact solutions are found which 
represent for instance periodic and quasi-periodic solutions and standing and travel- 
ing waves. Some of these results have been found previously by other authors via a 
study of mode coupling coefficients and recently also by investigating 'bushes of normal 
modes'. The method of this paper is similar to the latter method and much more sys- 
tematic than the former. We arrive at previously unknown results without any difficult 
computations. It is shown moreover that similar invariant manifolds exist also in the 
Klein-Gordon lattice and in the thermodynamic and continuum limits. 

1 Introduction 

The Fermi-Pasta-Ulam (FPU) lattice is a discrete model for a continuous nonlinear 
string, introduced by E. Fermi, J. Pasta and S. Ulam 0]. This string is modeled by 
a finite number of point masses which represent the material elements of the string. 
Each of the point masses is an oscillator that interacts with its nearest neighbors only. 
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Assume that the lattice consists of a finite number n G N particles. Define qj G M the 
vertical position of the j-th particle. We distinguish two different types of boundary 
conditions. We speak of fixed boundary conditions if the first and the last particle 
do not move, meaning that we have qo = Qn = Q for all time. The FPU lattice with 
fixed boundary conditions models a string with Dirichlet boundary conditions. It is 
also possible to choose periodic boundary conditions, in which case the first and the 
last particle are identified, that is qo = qn for all time. The FPU lattice with periodic 
boundary conditions models a circular string. Both types of boundary conditions occur 
very often in the literature. In this paper we shall only consider lattices with periodic 
boundary conditions, as it will become clear that each lattice with fixed boundary 
conditions is naturally embedded as an invariant manifold of an appropriate periodic 
lattice. The particles of the periodic lattice are labeled by elements of the cyclic group 
^/nZ- The Hamiltonian equations of motion for the FPU lattice are derived as follows. 

The space of positions q = (gi, . . . of the particles in the lattice is M". The 
space of positions and conjugate momenta is the cotangent bundle r*M" of M", the 
elements of which are denoted {q,p) = {qi, . . . ,qn,Pi, ■ ■ ■ ,Pn)- T*M.'^ is a symplectic 
manifold, endowed with the symplectic form dq A dp = Yl^=i ^^Ij ^ ^Pj- smooth 
function H : T*W^ M now induces the Hamiltonian vector field Xh given by the 
defining relation (dq/\dp){XH, ■) = dH. In other words, we have the system of ordinary 
differential equations qj = pj = — f^. 

The Hamiltonian function for the FPU lattice with periodic boundary conditions 
and n particles consists of a kinetic energy and a potential energy. The potential 
energy is assumed to depend only on the vertical distance between pairs of neighboring 
particles. Hence the Hamiltonian is 

^= E + ' (1-1) 

in which : M — R is a Lennard-Jones potential energy density function of the form 

W{x) = ^x' + ^x' + ^x'' . (1.2) 

The a,/? are real parameters measuring the nonlinearity in the forces between the 
particles in the lattice. We also write 

iez/„z ^ ^ 

In which 

^-(^) = h.^ ^"^^^ ~ 

is a polynomial in q of degree m. 

For a = /? = 0, the Hamiltonian is quadratic and the equations of motion linear. 
The solutions can be written down explicitly and the motion is in fact completely 
integrable. But for a,/? 7^ 0, the system is much harder to analyse. It can be inter- 
preted as a nonlinear perturbation of the integrable linear Hamiltonian system. Fermi, 
Pasta and Ulam expected that a many particle system such as the FPU lattice would 
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be ergodic due to these nonlinearities, meaning that almost ah orbits densely fill up 
an energy level set of the Hamiltonian H. Ergodicity would in theory have to lead 
to 'thermalisation', i.e. equipartition of energy between the various Fourier modes of 
the system. FPU's nowadays famous numerical experiment was intended to investi- 
gate how and at what time-scale this thermalisation would take place. The result was 
astonishing: it turned out that there was no sign of thermalisation at all. Putting 
initially all the energy in one Fourier mode, they observed that this energy was shared 
by only a few other modes, the remaining modes were hardly excited. Within a rather 
short time the system returned close to its initial state and thus behaved more or less 
quasi-periodically. 

The observations of Fermi, Pasta and Ulam greatly stimulated work on nonlinear 
dynamical systems. Nowadays people tend to explain the FPU experiment in two 
ways. In 1965 Zabuski and Kruskal ||l^ considered the Korteweg-de Vries equation as 
a continuum limit of the FPU lattice and numerically found the first indications for the 
stable behaviour of solitary waves. We now know that the Korteweg-de Vries equation 
is integrable. This clearly suggests an explanation for FPU's observations, although 
our understanding of the relation between the FPU lattice and its infinite dimensional 
limits has until now been quite disappointing. 

Another, possibly correct explanation for the quasi-periodic behaviour of the FPU 
system, is based on the Kolmogorov-Arnol'd-Moser (KAM) theorem. As is well-known 
||2|, the solutions of an n degree of freedom Liouville integrable Hamiltonian system 
are constrained to move on n-dimensional tori and are not at all ergodic but peri- 
odic and quasi-periodic. The KAM theorem states that most invariant tori of such an 
integrable system persist under small Hamiltonian perturbations, if the unperturbed 
integrable system satisfies a certain nondegeneracy condition which is called the Kol- 
mogorov condition. This shows that quasi-periodic behaviour can also be observed in 
certain nonintegrable Hamiltonian systems. Although several authors, starting with 
Izrailev and Chirikov Q, have therefore stated that the KAM theorem explains the 
observations of the FPU experiment, it has for a long time been completely unclear 
how the FPU system can be viewed as a perturbation of a nondegenerate integrable 
system. At least the linear Hamiltonian system for the FPU lattice with a = /3 = 
does certainly not satisfy the Kolmogorov condition. This gap in the theory was re- 
cently mentioned again in the review article of Ford and the book of Weissert 

The only rigorous results in this direction that are known to me were obtained by 
Nishida |12| and Rink |14]. These authors compute a Birkhoff normal form for the 
FPU lattice. This Birkhoff normal form is an approximation of the Hamiltonian of the 
FPU lattice. Although the results in |12] are unfortunately incomplete, it is proven in 
1 14] that this Birkhoff normal form is indeed often nondegenerately integrable. This 
explains why the FPU lattice at low energy can not be ergodic but must display a lot 
of quasi-periodic behaviour. On the other hand, many numerical studies indicate that 
above a certain energy threshold the lattice indeed thermalises. Reference [^] contains 
a rather complete overview of these results. 

Contrary to these more or less global approximation results, several authors have 
been trying to find exact low dimensional invariant manifolds for the FPU lattice. First 
of all because they represent interesting classes of solutions such as periodic and quasi- 
periodic solutions and standing and traveling waves. But also because it is believed 
by some authors, see for instance Q, that the destabilisation of invariant manifolds 
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can lead to chaos and hence maybe to ergodicity. The present paper is inspired by 
the idea that the results of which were primarily obtained for FPU lattices with 
periodic boundary conditions, will to a large extent also be applicable to subsystems 
of these periodic lattices. We will see for instance that every FPU lattice with fixed 
boundary conditions can be viewed as such a subsystem. This paper does not contain 
any explicit KAM statements though, but focusses on finding invariant manifolds only. 

Most of the invariant manifolds that are known in the FPU lattice were discovered 
more or less emperically. In their original paper Fermi, Pasta and Ulam already 
remarked that if the nonlinearity coefficient a in (1^) vanishes and initially only waves 
with an odd wave number are excited, then waves with an even wave number will 
never gain energy. Later on, other invariant manifolds were discovered by studying 
mode coupling coefficients in detail, see for instance and In these papers it is 
shown that certain sets of normal modes will not be excited if they initially have no 
energy. 

As will be explained in section ^ studying mode coupling coefficients can be quite 
unsatisfactory. A more systematic method for finding invariant manifolds in a physical 
system should be based on the symmetries of this system. The only reference that 
exploits these symmetries for the FPU lattice is |5| in which so-called 'bushes of nor- 
mal modes' are computed. These 'bushes' are simply invariant manifolds of a certain 
type. Their definition and how to find them are discussed more elaborately in . The 
basic idea is the well-known physical principle that the fixed point set of a symmetry 
forms an invariant manifold for the equations of motion. In [^, several previously un- 
known 'bushes' are classified. After computing the irreducible representations of the 
symmetry group of the FPU lattice and introducing appropriate 'symmetry-adapted 
coordinates', the computation of these 'bushes of normal modes' is fairly simple. 

In the present paper it will be shown that the previously mentioned invariant mani- 
folds and many others can be found even without introducing Fourier modes and study- 
ing mode coupling coefficients and without computing irreducible representations and 
symmetry-adapted coordinates. We only have to compute the fixed point sets of the 
various symmetries. As we incorporate more symmetries then H] in our considerations, 
we find various invariant manifolds that were not discussed before, in particular for the 
so-called /3-lattice. Moreover, our results are not only valid for the FPU lattice, but 
for any lattice with the same symmetries, such as the Klein-Gordon lattice |11]. They 
also apply in the thermodynamic limit as the number of particles grows large, and in 
the continuum limit: we can point out several infinite dimensional invariant manifolds 
for a rather broad class of nonlinear homogeneous partial differential equations. Some 
of them have been found previously in a very unpractical way. 



2 Quasi-particles 

Since we want to be able to compare our results with previous work, we introduce 
Fourier modes in this section. These Fourier modes are at the same time the 'symmetry- 
adapted coordinates' of [^. It is natural to view the solutions of the FPU lattice as a 
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superposition of waves and to make the following Fourier transformation: 

^i = ^ Yl ^^<ik (2-1) 



Pi 

/n 



1 r — ^ 'Z-Kijk 

-= ^~^P^ (2-2) 



1 2;rijfc j 1 if j = mod n 

n ' 1 if j 7^ mod n 



Using that 



one easily calculates that {qj-,qk} = {Pj-,Pk} = and {qj-,Pk} = ^jk, the Kronecker 
delta. Hence, are canonical coordinates. They are traditionally called phonons 



or quasi-particles. Written out in phonons, the FPU Hamiltonian (LI) reads as follows. 
The kinetic energy becomes: 

Yl = \Pn + \P\ + Yl PlP^-i ' 

where it is understood that the term occurs only if n is even. The potential 

^ 2 

energies Hm (for m = 2,3, 4) become 



-y (- y 

,!n 2 \ / ^ , 



27ri(j + l)fc 27rijfc , 

m! I v/n ^ 

jez/„z \ fcez/„z 



S^fct*^— Omod n 



in which the sum is taken over multi-indices S Z" for which \6\ := \6k\ = m. We 
also used the multinomial coefficient (™) := nfl^- have obtained a rather compact 
and tractible formula for the Hamiltonian in phonon-coordinates. 
Let us also introduce real-valued phonons. For 1 < k < ^ define 

Qk^ {qk + Qn-k)/ V2 = y - y^cos(^^^)gj , Qn-k = i{qk~ qn-k)/V2 = y - y^sin(^^^)gj 

Pk^{Pk+Pn-k)/V2 = W- V'cos(^:^)pj , Pn-k=i{pn-k-pk)/V2 = -\/ - sin(^:^)pj 
V n -"^ — ' n \ n ^ — ' n 
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and 

^ n - n 



Qn ^ 9n = 9j 1 Pn ^ Pn = 



The transformation {q,p) ^ (Q, P) is again symplectic and one can express the Hamil- 
tonian in terms of Q and P. In the case that a = /? = 0, that is for the harmonic FPU 
lattice, one gets 

iez/„z i=i 

in which for j = 1, . . . , n the numbers coj are the well-known normal mode frequencies 
of the periodic FPU lattice: 

jvr 

UJj := 2 sin{ — ) 
n 

Note that written down in real- valued phonon coordinates, the equations of motion of 
the harmonic lattice are simply the equations for n — 1 uncoupled harmonic oscillators 
and one free particle. The situation is not so simple anymore li a, (3 ^ 0, when the 
normal modes interact in a complicated manner that is governed by the Hamiltonians 



n-l 



e:\e\=m k=l 



in which the cq are certain coefficients. An expression for the cg can in principle be 
obtained from the formulas for the Hamiltonian Hm{q) and the mapping q ^ Q. For 



instance Hi{Q) can explicitly be found in [13|, although its computation is not given 
there. 

Note that H is independent of Qn = qn = Ylj 1j- Hence the total momentum 
Pn = Pn = Y^j Pj is a constant of motion and the equations for the remaining 
variables are completely independent of {Qn,Pn) = iPn,qn)- It is common to set the 
latter coordinates equal to zero, or to neglect them completely. Thus one removes the 
total momentum from the equations of motion. Equivalently, one could also perform 
the Marsden-Weinstein reduction of the symmetry induced by the flow of Xp^ = ^|-, 
cf. g or 0. 



In the nonlinear system, the other normal modes interact in a complicated manner, 
even though not every possible coupling term occurs. Only those monomials q^ = 
rifc^fc* present in Hm{q) for which "Ylik^^k = modn, whereas H^iQ) contains 
only the monomials = Q^j!' for which cq ^ 0. In the next section we will see 
that this is a consequence of discrete symmetries in the system. 

It is exactly the fact that not every coupling term occurs which accounts for the 
existence of various invariant manifolds, see § and m. Let A Cl/ni- Then the 
manifold spanned by modes in A is 

:= {(Q, P) e r*M"|Q,- = Pj = Q Vj ^ A} . 
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In several cases, these M;f are invariant manifolds for the equations of motion. In Q 
and they are then called 'bushes of normal modes'. We will not use this terminology. 
One readily infers from the equations of motion Qj = , Pj = — -^q- that M;^ is an 
invariant manifold (a 'bush') if and only if =0 for all 6 with the property that 9j = 1 
for some j ^ A and 6k = for all k ^ Au{j}. Making use of this fact, several invariant 
manifolds have been discovered. If n is even, one can for instance choose ^ = {^}. It 
is then obvious that A satisfied the required property since j + {m — 1)^ ^ mod n. 

{-} 

The solutions in the invariant manifold M„ ^ are of the form qj{t) = ^ (t)- This 
type of periodic solutions in which neighbouring particles are exactly out of phase, is 



well-known. In [13| a linear stability analysis is given for this solution in the /3-lattice 
(i.e. a = 0) and in |^ a similar linear stability analysis is given for this solution in the 
a-lattice (i.e. /3 = 0). 

Studying mode coupling coefficients in this way, several invariant manifolds have 
been discovered. In |Q it is shown that if a = and n is even, Mi''"'-'"^ and 
MA ' '•••'""^ are invariant. Poggi and Ruffo [||] show that M^^' a ^nd M^^ - 4 ^ 
are invariant. 

The above method is rather simple and easily understood but has the following 
limitations: 

1. An explicit expression for the cg is required. 

2. The method becomes more elaborate if one wants to find invariant manifolds of 
higher dimensions. 

3. There is no a priori 'physical' reason why a certain M:^ will be invariant. 

4. Invariant manifolds might exist that are not of the form M;f for some A C 'Z/nZ- 

5. It is not clear whether the discovered invariant manifolds will also be present in 
the continuum limit or in other one-dimensional lattice systems. 

For these reasons, studying mode coupling coefficients is rather unsatisfactory. With 
the method presented in the following sections of this paper it is possible to detect 
easily many more invariant manifolds. They arise in a natural way as fixed point sets 
of symmetries. 



3 Symmetry 

The Hamiltonian function (1.1) of the periodic FPU lattice has discrete symmetries 



with important dynamical consequences. Let us discuss symmetries in general here. 
Assume that P : T*R" — > T*R" is a linear isomorphism with the following two prop- 
erties: 

1. P is symplectic, i.e. P*{dq A dp) = dq A dp. 

2. P leaves the Hamiltonian invariant, i.e. P*H := H o P = H. 

Under these assumptions, P is called a symmetry of H. The set of symmetries of H is 
a group under composition of functions. This group is denoted Gh- 

For every symmetry P £ G we find that the Hamiltonian vector field Xh induced 
by H is equi variant under P: P*Xh = Xp*H = Xh- In other words: if 7 : M — > T*W^ 
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is an integral curve of Xh, then P o 7 : M — s- T*W^ is also an integral curve of Xff. 
This implies that P commutes with the flow of Xh , that is e'''^^ o P = P o e*^^^ . 
Of particular dynamical interest is the fixed point set of a symmetry P, 

Fix P = {{q,p) e T*R-\P{q,p) = (3.1) 

Let {q,p) e Fix P, then P{e^^"{q,p)) = e^^"{P{q,p)) = e^^"{q,p). So Fix P is 
an invariant manifold for the flow of Xh- This explains why fixed point sets are so 
interesting. When G' C Gh is a subgroup, then a fixed point set is also defined for 
it: Fix G' = ripgc/Fix P. These are of course also invariant manifolds and they are 
commonly studied. 

Fixed point sets of symmetries and fixed point sets of subgroups have a very simple 
relation. When Pi, . . . , Pm are symmetries, then (Pi, . . . ,Pm) C Gh is the smallest 
subgroup of Gh containing Pi, ... , P^. The symmetries Pi, ... , Pm are called gener- 
ators for this subgroup. One readily checks now that Fix (Pi, . . . ,Pm) = n^Fix Pj. 
It therefore suffices to study the fixed point sets of seperate symmetries. If G' is a 
subgroup of G that is generated by the symmetries Pi, ... , P^,, then the fixed point 
set of G' is simply the intersection of the fixed point sets of the seperate symmetries 
Pi, ... , Pm- As the number of elements of Gh can be considerably less then the number 
of subgroups of Gh-, we prefer to study fixed point sets of seperate symmetries first 
and take their intersections later. 

Let us discuss the symmetries of the FPU lattice now. Define the linear mappings 



P, S, r : ^ W 


by 










R ■- 


(91 : 


,Q2,-- 


■ 5 Qn-l,Qn) 1- 


(92,93, • • • ,9n,9l) 


(3.2) 


S : 


{Q1, 




• ) Qn-liQn) 1- 


(-9n-l, -9n-2, • • • , -91, -9n) 




T : 




■ q2,-- 


• 1 Qn-l-iln) 1- 


(-9l,-92,---,-9n-l,-9n) 





The mappings (g,p) ^ {Rq,Rp), {q,p) ^ {Sq,Sp) and {q,p) ^ {Tq,Tp) from r*]R" 
to r*M" are also denoted R, S and T respectively. They satisfy the multiplication 
relations R^ = S"^ = T"^ = Id and RS = SR~^, while T commutes with everything. 
Hence the discrete group (P, S) := {Id, P, P^ . . . , P"-\ S,RS, R^-'^S} is a rep- 
resentation of the n-th dihedral group the symmetry group of the n-gon, whereas 
(P, S, T) is a representation of Dn x 

R, S and T are symplectic maps and R and S leave the Hamiltonian H invariant. T 
only leaves H invariant if the potential energy density function W is an even function, 
in other words if a = 0. When a ^ 0, then Gh = (P, >S') is the symmetry group of H, 
whereas Gh = (P, -S", T) if W is even. 

In the coming sections we shall investigate the various invariant manifolds Fix P for 
P G Gh- We shall describe them in terms of the original coordinates {q,p), but also 
in phonon-coordinates {q,p) and {Q,P)- Therefore it is interesting to write down how 
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R, S and T act in complex phonon coordinates: 



R:{qi, 


92, •• 


• ,9n- 


-1, 


9n) 






.P2, • • 


' • ,Pn- 




>Pn) ^ 




S -{qi, 


92, •• 


■An- 


-1, 


9n) 


( — 9n-l, —qn-2, • • • , —91, — 9n) , 




.P2, • • 


■ • ,Pn- 




>Pn) ^ 


(-Pn-1, -Pn-2, • • • , -Pl, -Pn) ■ 


T -.{qi, 


92, •• 


■An- 


-1, 


9n) 


{-qi, -q2, ■ ■ ■ , -qn-1, -qn) , 




,P2, • • 


■ ■ ,Pn- 




,Pn) ^ 


{-pl, -p2, -Pn-1, -Pn) ■ 



Note that by performing the transformation to complex phonons, R has been diago- 
nalised, whereas the actions of S and T have not at all changed. This means that {q,p) 
are what in Q is called 'symmetry-adapted coordinates'. They are actually adapted 
to the subgroup (R) of Gh- The action of on a monomial is also very simple: 

R* (Hq'A =e2-S.fce./n Ylql- 
\ k / k 

In other words, the monomial is i?-symmetric if and only if kOj. = Omodn. So 
i?-symmetry is the reason why only these monomials occur in the FPU Hamiltonian. 



4 Invariant manifolds for arbitrary potentials 

In this section we study the invariant manifolds that are formed by the fixed point sets 
of elements of Gh = {R, S) = Dn- So it is not yet assumed that the potential energy 
density function W is even. 

For integers n and k, let gcd(n, k) be the greatest common divisor of n and k. For 

kez, 

Fix R'' = {qj = qj+gcd{n,k),Pj = Pj+gcd(n,k) Vj} 

is an invariant gcd(n, A;) degree of freedom symplectic submanifold of r*M". The 
Hamiltonian function -f^lpj^ j^k on the symplectic submanifold Fix R^ obviously simply 
models the periodic FPU lattice with gcd(n, k) particles. In this way, the periodic 
lattice with k particles is naturally embedded in the lattice with n particles if k divides 
n. In phonon coordinates. 

Fix R'' ={q, =Pj=0 Vi ^ mod 

gcd(n, k) 

71 

So if k divides n, then Fix R = Mn^' ^ ' "' * ' and is hence spanned by modes 
which represent a repeating spatial pattern with period k. 

{ — n} 

If for instance n is even, then Fix R? = Mn ^ ' is the two degree of freedom invari- 
ant manifold spanned by the §-th and the n-th normal modes. If we as usual neglect 
the n-th mode, which moves independently of all other modes, we find that Fix R? 
consists of all solutions of the form qj{t) = Q!i{t). These are the previously men- 
tioned periodic solutions in which neighboring particles are exactly out of phase. On 
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the other hand one has for even n that Fix R2 = M„ ' ' ■ ' . It consists of all even 
modes. ^ 

If 3 divides n, then Fix = Mi^'^'"\ whereas Fix = jv/i^'^'-'""^'"^ 
Etcetera. These invariant manifolds were discussed already extensively in 

The following invariant manifolds are only briefly discussed in For arbitrary / G Z 
we can study 

Fix R'S = {qj = -qi-j,Pj = -pi-j Vj} = {qj = -e ^Qn-j , Pj = -e^pn-j Vj} = 

llTT llTT llTT llTT ft 

{Qj cos(^) + Q„_j sin(^) = cos(^) + P„_, sin(^) = V 1 < j < - , 
n n n n 2 

It is a (2n — 2 — (—1)' — (— 1)"^')/4 degree of freedom symplectic subspace of r*M". 
Note that Fix R^S is not always of the form M;^ for some A. On the other 

hand, Fix S = mJ ^^^^"'^ and if n is even, then Fix R2S = Mn^'"' ^'^'^ ^'"'^ = 
^{i|2_j_2. i imod 2}u{i| 2 <j<'T., i Omod 2}^ ^^^^ instance for n = 8 these are Mg^'^''^^ 

and Mg^'-''^^ 

If both n and / are even, then Fix R'^S has dimension n/2 — 1 and in Fix ii'S* we 
have qi = qn+i = 0. In other words, if n is even, then for every even / the Hamiltonian 

2 2 

function H\p^^ on the symplectic subspace Fix R'^S models the FPU lattice with 
fixed boundary conditions and n/2 — 1 moving particles. Hence, the FPU lattice with 
fixed boundary conditions and n/2 — 1 moving particles is naturally embedded in the 
periodic FPU lattice with n particles. This is the reason why we do not study FPU 
lattices with fixed boundary conditions separately. 



5 Invariant manifolds for even potentials 

If the potential energy density function W is even, then also T is a symmetry and the 
full symmetry group of the FPU Hamiltonian is Gh = {R,S,T) = Dn x Z/2z- Let us 
study the fixed point sets of the symmetries R^T and RI^ST which have not yet been 
discussed in the previous section. Most results in this section are new, as the symmetry 
T was not considered in ||5|. 
For A; G Z, 

Fix R^T = {qj = -qj+gcd{n,k):Pj = -Pj+gcd{n,k) ^j} 

which is nontrivial only if n/ gcd(n, k) is even -and hence n must be even. In this case 
it is a gcd(n, k) degree of freedom invariant symplectic manifold. In phonons, 

Fi, fl'T =(5, = B = V, ^ ^^j^ mod ^j^} 

71 71 

= {Q . =R. = Vi / — — mod — - — -} . 

^ ^ 2gcd{n,k) gcd(n,/c)^ 

r 71 3n (2fc-l)n 

So if 2A; divides n, then Fix i?'=r = Mn^fe'^fe'-' 2fc \ 

The special choice k = ^ gives us the invariant manifold Fix R2T = Mn^'^'^''"'^ 
of all odd normal modes that was already discovered by Fermi, Pasta and Ulam 0. 
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The choice k = 1 gives us Fix RT = Mn^, the well known §-th mode. 

If n is divisible by 4, then Fix RiT = ^f^^'^'^^ ig invariant. This is a new 

^ n_ 3n \ 

result. The invariant manifold Fix R^T = Mn'^' * is discussed in [^]. It contains 
quasi-periodic solutions. 

For an n divisible by 6 we find the invariant manifolds Mn^'^'^^''"'"" and Mn " ' ^ ' ^ \ 
Etcetera. 

For / G Z, 

Fix RST = {qj = qi^j,Pj = pi^j Vj} = {Qj = e " Qn-j , Pj ^ e " Pn-j Vj} = 

{Q . siii(^) - Q„_, cos( — ) = Pj sin(^) - cos(^) = V 1 < j < - , 
n n n n I 

=(-l)'Q. , p. =(-l)'P.} 

is an (2n — 2 + (—1)' + (— 1)"^')/4 degree of freedom invariant symplectic manifold. 
Note that again Fix R^ST is not always of the form M;^, but that on the other 

hand Fix ST = Afi'''°-'-^^ and if n is even, Fix R^ ST = Mi°'"-''2'"-=^'^'-> = 

,.{i|0<j<f , i=Omod 2}U{i|f<j<n, j=lmod 2} or-- r r, ,.{1,2,3,4} 

Mn .bo tor mstance tor n = 8 these are Mg 

and Mt''''''^. 



6 Examples of intersections 

We have studied the fixed point sets of the elements of the symmetry groups Gh = 
(i?, S) and Gh = {R, S, T). They are equal to the fixed point sets of subgroups of Gh 
that are generated by one element. A fixed point set of a subgroup generated by more 
than one element must be the intersection of some of the fixed point sets that were 
already discussed. We will give just a few examples here. 



If 3 divides n, then Fix R^ n Fix S" = M„ ^ \ whereas Fix R^ n Fix ST 
The latter is only invariant if the potential W is even. 



r 3n \ / H H \ 

If 4 divides n, then Fix fi^ f^ S = M}, " \ Fix i?^ n Fix ST = ' ^ ^ 



and 



Fix i?27^n Fix ST = Mn" ■ 

If 5 divides n, then Fix i^^nFix S = " ' " \ whereas Fix ii^pFix ST = M^''' \ 
If 6 divides n, then Fix i?'' n Fix S" = ^ ' e ^nd Fix i?^ p Fix ST = « ' ^ ' ^ \ 

And we find that Fix i23r = Mn^'^' « ^ can be split into Fix R^T n Fix S" = M„ « ^ 

and Fix R'^T n Fix ST = M„ " ' ^ . The normal mode solutions for the /^-lattice that 
{— } 

lie in M„ ® have as far as I know never been discussed in the literature. 

One can proceed and compute, if k divides n, the intersections of the various fixed 
point sets of R , S, R2 S, R T, ST and R2ST. We choose not to make a systematic 
classification of the results, since most invariant manifolds in the FPU lattice are not 
even of the form M:^ for some A. 
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7 Other lattices and the continuum hmit 



A major advantage of our method is that fixed point sets of symmetries are invariant 
manifolds in any Hamiltonian system admitting these symmetries. Hence we expect to 
find the invariant manifolds that we discovered in the FPU lattice with periodic bound- 
ary conditions also in other one-dimensional spatially homogeneous lattices, such as 
the Klein-Gordon (KG) lattice [11|. The KG lattice with periodic boundary conditions 
has the Hamiltonian 

in which is a potential energy density function. The KG lattice models a one dimen- 
sional mono-atomic structure with small coupling between the atoms. It is clear that 
the mappings R and ST, see formulas (3.2), again leave this Hamiltonian invariant, 
whereas S and T separately have this property if W is an even function. Thus 
we have again found symmetries and their fixed point sets are invariant manifolds. In 
particular, the invariant manifolds that we discovered in the FPU lattice with even 
potential are also present in the KG lattice with even potential. 



In the thermodynamic limit it is assumed that the FPU lattice consists of a count- 
ably infinite number of particles, labeled by j € Z. The equations of motion are 
Hamiltonian equations on T*M^ with Hamiltonian 



2 

The symmetries are now induced by 



R 
S 
T 



{. . . ,q-i;qo,qi,. . .) ^ (. . . , go; 9i, • • •) 
{...,q-i;qo,qi,...) ^ {. . . , -qi, -qo; -q^i, . . .) 
{...,q^i;qo,qi,...) (. . . , -go, -ft, • • •) 



The finite dimensional manifold Fix i?" models an infinite lattice with a spatially re- 
peating pattern of period n. Or, equivalently, the periodic lattice with n particles. 
Inside Fix i?" we find again the invariant structures that were discussed previously in 
this paper. The invariant manifold Fix R'^S is an infinite dimensional one. It consists 
of solutions with qj = —q-n-j that are anti-symmetric around j = n/2. Etcetera. Sim- 
ilar conclusions hold of course for the thermodynamic limit of the KG lattice. 



Our results are also valid in the continuum limit, when the discrete lattice equations 
are replaced by a homogeneous partial differential equation. Consider for example for 
X G M/z the equation 

utt = u^x + f{u) , 

for / : M ^ M. This equation can also be written as the system of equations 

ut = v, vt = Uxx + f{u) , 

which have the Hamiltonian 

H= j \v{xf + ]-u^{xf - F{u{x)) dx 
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in which F' = f. Define the symplectic operators 



7^" : u{-) ^ u{a + ■), v{-) ^ v{a + ■) 
S : u{-) ^ —u{— •), f(-) I— > —v{— •) 
T : u{-) I— > — u(-), f (•) I— > — v(-) . 

The constant a G M/z is arbitrary. Clearly, H is invariant under TZ"^ and 5T. is 
invariant under TZ"", S and T separately if and only if F is even, that is if and only if 
/ is odd. 

The fixed point sets of these symmetries are invariant manifolds, possibly of infinite 

dimension. If a ^ Q, then Fix TZ"" consists of constant solutions only, but if a = - 

p 

is rational and gcd{p,q) = 1, then Fix TZi represents the solutions with u(t,x) = 

1 - 1 

u{t,x+-). Fix 7^9 T consists of solutions with M (re) = — The latter is nontrivial 

only if q is even. For arbitrary a, Fix TZ^S contains solutions with u{x) = —u{a — x) 

and Fix IZ'^ST represents solutions with u{x) = u{a — x). 

It is natural to use the Fourier transformation 

u{x,t) =Y^Uk{t)e'^^^ , v{x,t) = Y^Vk{t)e'^-^ 

and to express the fixed point sets in terms of the Fourier variables {uk,Vk)kGZ- We 
then find for instance the following invariant manifolds 

Fix TZ^ ={uk = Vk = Oyk^O mod q} = m^- -9,o,g,2g,...} 
Fix TZ^T = {uk = Vk = Oyk^q mod 2q} = m^- -m.Sq,...} 

Etcetera. 



[p^], [15 1 and |jT^ study the equation uu = u^x + by the Galerkin-averaging 
method. By an analysis of mode coupling coefficients they discover that the manifolds 
j^{...,-2(jr,-g,o,g,2(3,...} y^{...~Zq~q,q,Zq,..?i invariant in a certain finite dimensional 

system of differential equations, the Galerkin-averaging approximation, which approx- 
imates the original partial differential equation. We arrive here at the much stronger 
result that their conclusions hold for any odd nonlinearity / and in the original partial 
differential equation. 



8 Discussion 

In a systematic way we found various invariant manifolds for the FPU oscillator lat- 
tice with periodic boundary conditions. These invariant manifolds represent interesting 
classes of solutions such as periodic and quasi-periodic solutions, standing and traveling 
waves and embedded lower dimensional FPU lattices with periodic or fixed boundary 
conditions. They are moreover interesting since it is believed by some authors |Q] that 
destabilisation of these invariant manifolds can lead to chaos. Some of the invariant 
structures that we found have previously been discovered by other authors by an analy- 
sis of mode coupling coefficients. Our method on the contrary is similar to the method 
of 'bushes of normal modes' and looks for fixed point sets of symmetries which are 
natural invariant manifolds. We can derive our results without computing mode cou- 
pling coefficients explicitly. In fact, it is not even necessary to introduce normal modes 
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at all as an expression for the invariant manifolds can simply be obtained in terms of 
the original physical variables, which are the positions and momenta of the particles 
in the lattice. In this way, we find several previously undiscussed invariant manifolds 
in the FPU lattice. The same invariant manifolds are present in other homogeneous 
Hamiltonian lattices such as the Klein-Gordon lattice and even in lattices with an in- 
finity of particles. In the continuum limit, when the lattice equations are replaced by a 
homogeneous partial differential equation, we point out analogous infinite dimensional 
invariant structures. 
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